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An Interacting Dark Energy Model for the Expansion History of the Universe 
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We explore a model of interacting dark energy where the dark energy density is related by the 
holographic principle to the Hubble parameter, and the decay of the dark energy into matter occurs 
at a rate comparable to the current value of the Hubble parameter. We find this gives a good fit 
to the observational data supporting an accelerating Universe, and the model represents a possible 
alternative interpretation of the expansion history of the Universe. 
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I. INTRODUCTION 

Observations of type la supernova indicate that the 
Universe is making a transition from a decelerating phase 
to an accelerating one. The conventional explanation for 
this behavior is that a dark energy component is com- 
ing to dominate over and the matter-dominated phase is 
giving way to a phase dominated by a dark energy com- 
ponent. In fact a good fit is obtained from the observa- 
tional data by assuming the cosmological model (ACDM) 
involving a cosmological constant A and cold dark matter 
in about the ratios flm,o = 0.3 and f^A.o — 0.7. 

In this paper we show that one can describe the cos- 
mological data with a model of dark energy that in- 
cludes an interaction that effects a transition (decay) 
of dark energy into matter. The model incorporates 
a holographic principle^, 0] to determine the dark en- 
ergy density of the Universe. The principle relates 
the dark energy scale to the Hubble horizon which has 
been the subject of speculation for applying holographic 
ideas to cosmology^ iSl jfl, J., 3j ili JJlJ- This is a 
common choice for imposing holography on cosmology, 
and it has the most natural thermodynamic interpreta- 
tion. However, this conditon imposed on the dark en- 
ergy yields an equation of state that is matter- like [llj 
and therefore inconsistent with observational data. Sub- 
sequently efforts have been made to impose a holo- 
graphic constraint based on some other physical horizon 
such as the particle horizon 1 12l 131 or the future event 
horizon [T3.ll5[IT^IrRIT8l.h9l.l2Ctl2 j| or some other phys- 
ical condition (221 1231 bsl l26l l27r i2S 

A suitable evolution of the Universe is obtained when, 
in addition to the holographic dark energy, an interaction 
(decay of dark energy to matter) is assumed, and the de- 
cay rate should be set roughly equal to the present value 
of the Hubble parameter for a good fit to the expansion 
history of the Universe as determined by the supernova 
and cosmic microwave background (CMB) data. 

Interacting dark energy has been studied 
previously IM IIS III 111 111 Im III 111 111 111 111 E 113 



primarily with a goal of explaining the cosmic coinci- 
dence problem. A survey of the possible dynamics of 
dark energy can be found in Ref. |42l| . In the presence 
of an interaction the dark energy can achieve a stable 
equilibrium that differs from the usual de Sitter case 
or the approach to the stable equilibrium can be made 
more gradual. Such models offer the hope of solving the 
coincidence problem that exists in the ACDM model 
where there is no obvious reason why the transition from 
matter domination to domination by the dark energy is 
occuring during the current epoch. 

This paper presents a simple model that gives an ac- 
ceptable expansion history of the Universe in terms of 
a holographic condition on the dark energy that relates 
its size to the Hubble scale. This is a common choice 
for imposing holography on cosmology, and it has the 
most natural thermodynamic interpretation. The effec- 
tive equations of state of matter and dark energy coincide 
and behave like cold dark matter (CDM) at early times. 
The transition to behavior like a cosmological constant 
is effected by simply assuming there is a constant decay 
of dark energy into matter. The coincidence problem ap- 
pears in the interacting dark energy model as the choice 
of the scale for the dark energy interaction which must 
be close to the present value of the Hubble parameter. 



II. A FRAMEWORK FOR INTERACTING 
DARK ENERGY 



The continuity equations for dark energy and matter 



PA + 3i/(l -I- wa)pa 

+ 3H{1 + Wm)p m 



(1) 
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The interaction is given by the quantity Q, and the 
energy-momentum tensor remains conserved as long as 
the same factor Q appears on the right hand side of each 
equation. To appeal to our intuition and to facilitate 
comparison with the more familiar case of ACDM model, 
it is useful to define effective equations of state that play 
the role of the native equations of state when an interac- 
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tion is present. Following Ref. [20j, if we define 



r 

3H 



7^ • (2) 



where T — Q/pA is a rate, the continuity equations can 
be written in their standard form 



PA + 3H{l + wf)pA = , 
Pm + iH{\ + W''^)pra = 



Define as usual 



Sirpn 



If we restrict our attention to the flat case. 



one obtains a differential equation 



cIQa 
dx 



= -S^Ail - ^a) [wf 



..cffl 



(3) 



(4) 
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where x = ln(a/ao) and a is the Friedmann-Robertson- 
Walker (FRW) scale factor, and the "0" subscript here 
and on other quantities indicates the parameter value at 
the present time. For the ACDM model, where w'^ = 
WA = — 1 and w'^ = Wm — 0, the solution is the familiar 
one which can drive a small vacuum energy solution from 
a fixed point at J^a ~ toward a de Sitter Universe 
fixed point at Qa = 1. In the presence of a holographic 
principle and/or an interaction, the effective equations of 
state can depend on x, and there is the possibility that 
another fixed point can develop at a value of J^a other 
than or 1. 

As shown in Ref. (4(1] two physical assumptions are re- 
quired to determine the evolution of VIa and These 
can be chosen from the following list: (a) a "holographic 
principle" which specifics pA in terms of some length scale 
associated with a horizon, (b) an assumption for the in- 
teraction r, or (c) an assumption about the native equa- 
tions of state WA and Wm- The assumptions commonly 
used for (a) are assumptions for the length scale La in 



PA 



(7) 



Here the constant c represents an order one constant. 
Some choices have been to identify the length scale with 
a physical length such as the Hubble horizon, the parti- 
cle horizon, or the future event horizon, or perhaps some 
other parameter. The interaction rate in (b) can be speci- 
fied by an assumption of how the ratio of rates T/H varies 
as a function of Qa- The assumption for the equations of 
state in (c) have been largely confined to constant ones. 

It is important to note that specifying two of the three 
conditions determines the third. For example a generic 
choice for the length scale La and for the interaction rate 



F is inconsistent with a constant equation of state. It is 
easy to show that[40l| 



T = 3H{-1-wa) + 2 



La 
La 



(8) 



In addition to providing the connection between the var- 
ious rates and the equation of state, this also establishes 
that the interaction should be of the size of H. 

Using the above formulism one can show that (assum- 
ing Wm = 0) 



off 1,2 La 
= -'+3H-A- 

F flA 



(9) 

3H{i-nA) ■ ^^^^ 

Then one obtains from Eq. (j^Jl the differential equation 



wf = - 



dflA 

dx 



3r!A(l-^^A) 



J_La 
3H La 



3H (1 - flf 



(11) 



This evolution equation is general and is obviously most 
useful when an assumption for the holographic princi- 
ple and for the interaction are specified. For a constant 
equation of state together with a holographic principle 
the evolution equation becomes 



dflA 
dx 



^ SQa 



1 -I- WA^A 



_2_La 
2,H La 



(12) 



On the other hand if a constant equation of state and 
the interaction is specified, one can eliminate the explicit 
dependence on La/ La and obtain 



dilA 
dx 



3r!A(l-r!A) 



1 



-WA 



3H (1 - Qa) 



(13) 



by 



Finally the evolution of the Hubble parameter is given 



1P ~ H~dx ~ ^2^ 



Equations Ijllll and 114|l then determine a set of equations 
that give ^a{z) and L[{z) where z is the redshift, and 
include the usual ACDM model as a special case. 



III. INTERACTING DARK ENERGY MODEL 

The ACDM model is defined in terms of the preceding 
formalism as 



r = . 



(15) 



The equations of state of matter and dark energy are 
and —1 respectively, and the model predicts an evolution 
of the Hubble parameter as 



H{z)=Ho [{l+zfnrafi+nA,0 



1/2 



(16) 
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where flm,o and JIa.o are the present values of matter and 
dark energy. A good fit to the supernova data is given 
by flmfi ~ 0.3 and JIa.o ~ 0.7. 

Various physical assumptions have appeard for the 
holographic condition on L\/ L\. The cosmological con- 
stant (ACDM) requires L\/L\ = 0. Other conditions 
involve invoking a holographic argucmcnt such as relat- 
ing the length scale to various physical horizons. For 
example, one can relate La to the Hubble horizon (HH) , 
the particle horizon (PH) , or the future horizon (FH) for 
which 



^PH 
i?FH 



1/H 



* dt 



dt 



° da 



da 



(17) 



These yield 




The HH condition implies 

d^A 
dx 



= . 



(18) 



(19) 



so that (and is constant. This scenario has been 
shown to be inconsistent with the observational data|ll| 
in the absence of any interaction. However as we now 
proceed to show, if a constant rate of decay F of dark 
energy into matter is assumed, one can recover a good fit 
to the data. 

The PH and FH conditions have been considered more 
recently, and have been exploited^, |23| to attempt to 
achieve a fixed point solution for JIa- It has been shown 
that the FH condition is needed if one is to obtain a 
value of r^A that is suitably close to the valued obtained 
from fits to the supernova data. However assuming that 
f^A has acquired a value suitably close to its fixed point 
value does not guarantee a good fit to the data. In fact, 
by assuming PH or FH as the holographic condition and 
an interaction as a function on JIa, the quantity H /H^ 
would be a function of JIa alone as in Eq. (|14|) . If also 
Oa is approximately constant because it is near its fixed 
point, then the universe is characterized by two compo- 
nents with constnat equations of state and w^^ ■ 

Consider the case where a holographic condition asso- 
ciates the dark energy length scale with the Hubble pa- 
rameter (HH) and the interaction is set equal a constant 
of the same scale as the Hubble constant, 



La 



H 
H ' 



kHq 



(20) 



Here we require k to be an order one constant. As pointed 
out in Ref. the first condition in Eq. H2U|I in the ab- 
sence of any interaction (F = 0) gives rise to a descrip- 
tion of the universe inconsistent with the observational 
data. In this circumstance, the holographic condition 
implies that the cosmological evolution will appear as 
matter-dominated (w = 0) at all times. However when 
the interaction (decay) of the dark energy in Eqn. H2U|) 
is included, one can obtain a suitable fit. The Hubble 
parameter satisfies 

where tq — Qjn,o/^A,o- This evolution of the Hubble pa- 
rameter exhibits the needed features to agree with ob- 
servations: for large z a characteristic (1 -I- z)'^/^ be- 
havior of a matter-dominated era with an approach at 
smaller z to a constant H{z) characteristic of a de Sit- 
ter phase. A good fit is obtained when K/(3ro) w 0.62 
(see Section IV). Comparing to the ACDM solution in 
Eq. 1)16(1 , the best fits result for parameters roughly equal 
to riA.o = 2/3 and flm.o — 1/3 in the ACDM and simi- 
larly k/ (3ro) = 2/3 in the interacting dark energy model 
(the parameter k/(37'o) in the interacting dark energy 
model plays an analogous role to that of JIa.o in the 
ACDM). Nevertheless the best fit value for K/(3ro) in 
the interacting dark energy model is not equal to the 
best fit value of flA,o in the ACDM. 

The characteristics of the model can be understood 
from 



-1 



1 



dflA 



(22) 



H 3 

IP ^ 2 \ ^ ' SHl-flAj ' 2(1 -f^A) dx 

A constant Qa causes the last term to vanish and the 
l/r = riA/(l — ^^a) factor is constant and equal to its 
present value I/tq. Then 



H 



F 1 
3h7o 



-1 + — — 



(23) 



Then H{z) clearly exhibits a transition between a matter- 
dominated era {H / = —3/2) and a constant H era 
{H/H^ = when ^ Ki/o/(3ro) for T = kHq). At 
early times the interaction is negligible, and the holo- 
graphic condition enforces an expansion of the Universe 
that appears as matter-dominated. When the Hubble 
parameter becomes sufficiently small, the interaction be- 
comes effective and a fixed point solution is reached. 

While the interacting dark energy model predicts a 
constant Ha the coincidence problem emerges as the 
choice of scale F = kHq. The quantity k is an order 
one to give a good fit to all observational data. 

The holographic condition La/La = —H/H implies 
the effective equations of state are equal. 



„cff 



F 



3iJ (1 
F 



VtA) 



(24) 
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With the interaction F = kHq one recognizes that the 
common equation of state for the matter and dark energy 
components varies from to — 1 as the Hubble parameter 
decreases from large values in the early Universe to its 
asymptotic value {Hoo — kHq / (3ro)) . 



IV. QUANTITATIVE COMPARISON OF THE 
MODEL WITH OBSERVATIONAL DATA 

In this section we compare the interacting dark energy 
model predictions from Eq. H21() with supernova data. 
Comparison of the observational data with holographic 
models with noninteracting dark energy was performed 
in Refs. (21 and [l^ . The holographic condition based on 
the Hubble horizon (HH) is grossly inconsistent with the 
data as shown by Hsu[il|, but other holographic models 
based on the future event horizon (FH) can be made 
consistent with the data. 

The luminosity distance is defined as 



dL{z) = {l + z)H^^ 



dz' 



where E{z) — H{z)/Hq. In the interacting dark energy 
model, E{z) is obtained from the integration of Eq. 
The observational data for the supernovae is expressed in 
terms of an apparent magnitude and redshift. Assuming 
the supernovae have the same absolute magnitude M, 
then the extinction-corrected distance moduli is given by 



fi{z) = 51ogio(dL(z)/Mpc) + 25 



(26) 



Using the calculated distance moduli in Eq. for 
an interacting dark energy model and the supernovae 
data^^, one can perform a x-squared fit, 



E 



[A*obs(Zi) - f^mod{zi)Y 



(27) 
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z [redshift] 



FIG. 1: Comparison of the supernova data with the model 
prediction for fi{z) with ro/t = 1/2. The data shown in the 
Gold data from Ref . lial . The best fit occurs for ft=0.64 for 



(25) which the x is 179 for 157 degrees of freedom. 



of the interacting dark energy model is the same as the 
ACDM, the predicted value of shift parameter 



n 



1 



dz' 

W) 



(28) 



is easily accommmodated with a choice of tq/k sa 1/2 
and JIa^o ~ 0-7. The utility of the shift parameter TZ 
and the ^-parameter measured by SDSS is that they are 
independent of Ho 113 ■ Consequently the gross geometric 
features of this model agree with the expansion history 
predicted by the ACDM. An important question remains 
as to whether the detailed observations of the large scale 
structure can be made consistent with the matter density 
predicted at early times in this model (or some variant). 



where the sum runs over the supernova data points and 
where the ct^ are the experimental errors in each obser- 
vation. The ACDM model gives a of 178 for the 157 
data points in the Gold data sample li^ for 17^,0 = 0.27 
and SIa^o = 0.73. 

A comparison of the supernova data with the interact- 
ing dark energy model is shown in Fig. 1. Taking the 
Hubble constant as H^^ = 2997.9/1"^ Mpc, the best fit 
values for tq/k is 0.54 and for the Hubble constant is h 
=0.64. The Itr and 2a contours obtained when varying 
the parameters h and tq/k are shown in Fig. 2. 

The luminosity distance has dimensions of length and 
therefore scales with the inverse Hubble constant Hq. 
This means the model curve for ^(z) can be shifted up 
or down by adjusting Hq. The interacting dark energy 
model is also consistent with CMB measurement of the 
shift parameter '44', '45', and measurements of large 

scale structure by SDSS,4S|. Since the large z behavior 



V. SUMMARY AND CONCLUSIONS 

The ACDM model provides an explanation for all ob- 
servational data. However there remain a number of im- 
portant issues it must confront. Given our lack of un- 
derstanding of the dark energy, one can ask if there are 
other simple physical properties that dark energy might 
have that could equally well account for the data. 

In this paper we have shown that a model of dark en- 
ergy with a holographic condition relating the dark en- 
ergy to the Hubble parameter and a constant interaction 
of size roughly equal to the Hubble constant Hq can give 
a good fit to data. This model is characterized by a con- 
stant ratio of fim to JIa- Since the ratio of matter to 
dark energy is constant due to the holographic condition 
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view the evolution of the model as one comprised of one 
component whose effective equation of state varies be- 
tween when the decay of the dark energy into matter 
is negligible to an asymptotic value of —1 when the in- 
teraction become important. The transition between one 
regime and the other occurs when the interaction rate is 
comparable to the Hubble parameter. 

We have shown that the expansion history of the Uni- 
verse can be reproduced with a model of interacting 
dark matter. For the specific model presented here the 
functional dependence H{z) differs slightly between the 
ACDM and the model discussed in this paper. Future 
data from SNAP may be useful in discriminating between 
them. Finally information that goes beyond the recorded 
expansion history of the Universe may be useful for ruling 
out this kind of model. 



FIG. 2: The Ax^ = 2.30,6.17 contours {la, 2a) for the pa- 
rameters h and ro/n in the interacting dark energy model. 
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